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SUMMARY 


A three-dimensional implicit Full Navier-Stokes (FNS) analysis and a 3D Reduced 
Navier Stokes (RNS) initial value space marching solution technique has been applied to 
a class of separated flow problems within a diffusing S-duct configuration characterized 
as vortex-liftoff. Both Full Navier-Stokes and Reduced Navier-Stokes solution tech- 
niques were able to capture to overall flow physics of vortex lift-off, however more con- 
sideration must be given to the development of turbulence models for the prediction of the 
locations of separation and reattachment. This accounts for some of the discrepancies in 
the prediction of the relevant inlet distortion descriptors, particularly circumferential dis- 
tortion. The 3D RNS solution technique adequately described the topological structure 
of flow separation associated with vortex lift-off. 

NOMENCLATURE 


A t 

= 

inlet throat area 

C f 

= 

wall skin friction coefficient 

Dc 

= 

inlet throat diameter 

DC { * 

= 

distortion descriptor defined as the 

maximum ( Pt m — Pt^fq^ in any 60.0° sector 

L 

= 

length of inlet duct 

Me 

= 

inlet throat Mach number 

Pto 

= 

free stream total pressure 

Pt m 

= 

average total pressure at the engine face 

Pt m i„ 

= 

minimum total pressure at engine face 



in any sector of extent 60.0° 

^}ave 

= 

average dynamic pressure at the engine face ■ 

Re 

= 

inlet throat radius 

Ref 

= 

engine face radius 

Rey 

= 

Reynolds number based on throat diameter 

Tt o 

= 

free stream total temperature 

X , Y, Z 

= 

primary cartesian coordinates 

X* Yc, Zd 

= 

cartesian coordinates along inlet centerline 

6 

= 

boundary layer thickness 

P 

— 

fluid density 
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INTRODUCTION 


Modern tactical aircraft are required to be maneuverable at subsonic, transonic, 
and supersonic speeds, without giving up good cruise performance. Consequently ^proper 
integration of the engine inlet with the airframe is of paramount importance... Regarding 
the enhancement of inlet performance and operation, design for optimum airframe-inlet 
integration has the following goals: (1) to minimize approach flow angularity with respect 
to the inlet cowl lip, (2) to deliver uniform, high pressure recovery flow to the inlet face, 

(3) to prevent or minimize vortex, wake, and boundary layer ingestion by the inlet 
throughout the flight envelope, (4) to reduce FOD/hot gas ingestion by the inlet, and 
finally (5) to minimize the potential for flow field interference from weapon 
carriage/firing, landing gear deployment, tanks, pods, or other hardware. The combina- 
tion of inlet design and airframe integration must not only provide high pressure recovery 
to maintain the desired thrust levels, but also generate low flow distortion consistent with 
stable engine operation. 

Engine face flow distortion is one of the most troublesome and least understood 
problems for designers of modern inlet engine systems.'* 2 One issue is that there are nu- 
merous sources of flow field distortion that are ingested by the inlet or generated within 
the inlet duct itself. Among these sources are (1) flow separation at the cowl lip during 
maneuvering flight, (2) flow separation on the compression surfaces due to shock-wave 
boundary layer interactions, (3) spillage of the fuselage boundary layer into the inlet duct, 

(4) ingestion of aircraft vortices and wakes emanating from upstream disturbances, and 

(5) secondary flow and possibly flow separation within the inlet duct itself. Most aircraft 
have experienced one or more of these types of problems during development, particularly 
at high Mach numbers and/or extreme maneuver conditions, such that flow distortion at 
the engine face exceeded allowable surge limits. Such compatibility problems were en- 
countered in the early versions of the B70, the F-l 1 1, the F-14, the MIG-25, the Tornado 
and the Airbus A300 to name a few examples. 

Flow separation in diffusing S-duct geometries characterized as vortex lift-off has 
been studied by Harlof, Reichert, and Wellborn 2 using FNS analysis techniques and by 
Anderson and Farokhi 4 using an RNS approach. This class of separated flows are very 
common within inlet ducts, and are a major cause of pressure loss and distortion at the 
engine face, particularly circumferential distortion. Reduced Navier-Stokes solution 
technique using FLARE approximations 4 has been shown to physically describe the 
topological and topographical structure flow separation associated with vortex liftoff, 
however no detailed comparative study has been made between FNS and RNS viscous 
analyses for this phenomena over a wide range of inlet flow conditions. 

The overall objective of this study is to advance the understanding, the prediction, 
and the control of inlet distortion, and to study the basic interactions that are involved in 
the management of engine face distortion within inlet ducts using Computational Fluid 
Dynamics. This paper examines the phenomena of vortex lift-off using both Full 
Navier-Stokes and Reduced Navier-Stokes solution techniques, each with an algebric 
eddy viscosity turbulence model. Specifically, the goals of the present paper are: (1) to 
examine the capability of the Full Navier-Stokes and Reduced Navier Stokes analyses to 
describe the phenomena of vortex lift-off over a wide range of inlet flow conditions, (2) 
to characterize the phenomena of vortex lift-off and identify uncertainties in the analysis 
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of this interaction, and (3) to examine the prediction of the relevant inlet performance 
parameters of total pressure recovery and various engine face discriptors relative to ex- 
perimental measurements. 

THEORETICAL BACKGROUND 
Reduced Navier-Stokes Analysis 

Three dimensional viscous subsonic flows in complex inlet duct geometries are in- 
vestigated by a numerical procedure which allows solution by spatial forward marching 
integration, utilizing flow approximations from the velocity-decomposition approach of 
Briley and McDonald. 5 ' 6 The goal of this approach is to achieve a level of approximation 
that will yield accurate flow predictions, while reducing the labor below that needed to 
solve the full Navier Stokes equations. The governing equations for this approach have 
been given previously for orthogonal coordinates, and the approach has been applied 
successfully to problems whose geometries can be fitted conveniently with orthogonal co- 
ordinate systems. However, geometries encountered in typical subsonic inlet ducts cannot 
be treated easily using orthogonal coordinates, and this led to an extension of this ap- 
proach by Levy, Briley, and McDonald, 7 to treat ducted geometries with nonorthogonal 
coordinates. In generalizing the geometry formulation, Anderson, 8 extended the analysis 
to cover ducted geometries defined by an externally generated gridfile, such that it al- 
lowed for (1) reclustering the existing gridfile, (2) redefining the centerline space curve, 
and (3) altering the cross-sectional shape and area distribution without returning to the 
original gridfile. This version of the 3D RNS computer code is called RNS3D. The tur- 
bulence model used in RNS3D is that of McDonald and Camarata 9 which employes an 
eddy-viscosity formulation for the Reynolds stresses. 

Full Navier-Stokes Analysis 


The PARC3D 10 code, selected for this study, solves the full three-dimensional, 
Reynolds averaged, Navier-Stokes equations in strong conservation form with the Beam 
and Warming approximate factorization algorithm. The implicit scheme uses central 
differencing for a curvilinear set of coordinates. The code was originally developed as 
AIR3D by Pulliam and Steger. 11 Pulliam 12 later added the Jameson 13 artificial dissi- 
pation and called the code ARC3D. Copper 10 adapted the ARC3D computer code for 
internal propulsion applications and named the code PARC3D. The turbulence model 
used in PARC3D is the Baldwin- Lomax 14 model which is also an algebric eddy viscosity 
model. In the present computations, the turbulence model has been modified to improve 
the simulation of the reverse flow regions based on the study Deiwert. 15 In the regions 
of reverse flow, the inner layer is replaced with the with the outer model which extends 
all the way to the wall. In the absence of reverse flow, the conventional Baldwin-Lomax 
turbulence model is used. 
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RESULTS AND DISCUSSIONS 


The experiment on which this study is based comes from the AGARD Fluid Dy- 
namics Panel Working Group 13 numerical sub-group Test Case 3. 16 AGARD Working 
Group 13 was formed to investigate the subject of inlet performance using results from 
both experimental and computational techniques. Emphasis was placed on the assess- 
ment of inlet performance flow distortion, and the evaluation of relevant distortion pa- 
rameters over a wide range of inlet flow conditions. One of the test cases chosen by 
Working Group 13 was the RAE inlet model M2129, which is a circular inlet followed 
by an S-duct diffuser. The M2 129 inlet duct geometry, chosen by AGARD Working 
Group 13 because of the comprehensive data available, is shown in Fig. 1 and was based 
on a study by Willmer, Brown, and Goldsmith. 17 The centerline of the inlet defined in 
terms of the coordinate system shown in Fig. 1 is given by: 

Z c/ = -AZ^I-cos(*2k)] (1) 


where Xd is the x-coordinate of the inlet duct centerline, and A in the centerline offset. 
The radius distribution measured perpendicular to the duct centerline is given by: 



where R< is the inlet throat radius, K f )s the engine face radius, and L is the length of the 
inlet. For the purposes of the calculations, the M2 129 S-duct was nondimensionalized 
with respect to the throat radius, thus R, = 1.0, Kf= 1-183, L = 7.10, and A Z d = 2.13. 


A polar grid topology was chosen for the M2 129 S-duct which consisted of 49 ra- 
dial 49 circumferential, and 61 streamwise nodal points in the half plane, for a total 
number of 146,461 grid points. The internal grid was constnicted such that the transverse 
computational plane was perpendicular to the duct centerline. Grid clustering was used 
in the radial direction in order to redistribute the nodal points to resolve the high shear 
regions near the wall. The flow in the inlet was considered turbulent throughout. The 
inflow boundary layer condition corresponds to a shear layer thickness d/R, - 0.120, and 
were applied one diameter upstream of the inlet entrance in the constant area extension. 

A summary of the inlet test conditions used in this study is presented in Table I, 
and covers the inlet throat Mach number range from 0.1 to 0.794. It also includes the 
AGARD Test Case 3.1 an 3.2 test conditions defined by: 


Test Case 3.1 Test and Initial Conditions 


Total Pressure 
Total Temperature 
Throat Mach Number 
Throat Diameter 
Throat Area 

Reynolds Number (based on D t ) 


Pt 0 = 29.889 in. Hg 
Tt 0 = 293 ° K 
Mi = 0.794 
Di = 5.07 1 in. 

Ai = 25.254 in? 

Rey = 1 .848.x 10 6 
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Test Case 3.2 Test and Initial Conditions 


Total Pressure 
Total Temperature 
Throat Mach Number 
Throat Diameter 
Throat Area 

Reynolds Number (based on A,) 


Pt 0 = 29.865 in. Eg 
Tt 0 = 293° K 
Mi = 0.412 
A = 5.071 in. 

At = 25.254 in? 

Rey = 1 .1 58jc10 6 


For the purposes of examining the the separation characteristics within the M2 129 
inlet S-duct with both Full Navier-Stokes and Reduce Navier-Stokes solution techniques 
at a low inlet throat Mach number, a third Test Case was defines as 


Test Case 3.3 Test and Initial Conditions 


Total Pressure 
Total Temperature 
Throat Mach Number 
Throat Diameter 
Throat Area 

Reynolds Number (based on A) 


Pt 0 — 29.889 in. Eg 
Tt 0 = 293° K 

Mi = 0.200 

A, = 5.071 in. 

Ai = 25.254 in? 

Rey = 0.594*1 0 6 


Also presented in Table I is the computed location of separation X SfP from both the 
FNS and RNS analyses (as measured in terms of the arc length along the inlet centerline), 
the area averaged engine face total pressure recovery P R, and the AC® engine face dis- 
tortion, both determined from the flow values on the computational mesh. 

The initial flow Field for the Full Navier-Stokes computations was obtained from 
the Reduced Navier-Stokes (RNS) solution. This was done to reduced the computing 
time of the Full Navier-Stokes solver. In addition, the imposition of the inlet boundary 
layer, which was held Fixed during the time-marching FNS solution, was straight forward 
since* the initial flow Field had the correct initial inlet boundary layer. Examination of 
saving in the computing time of the Full Navier-Stokes computations using the RNS sol- 
ution as the initial flow Field compared to that using a simple uniform flow Field is cur- 
rently underway. 

Figures (1) and (2^ present the effect of inlet throat Mach on the area averaged 
total pressure recovery Pt e ffPto and DC® engine face distortion as determined from the 
Reduced Navier-Stokes and Full Navier-Stokes solution technique. The two solution 
techniques gave essentially the same total pressure recovery characteristics over the inlet 
throat Mach number ranged investigated, Fig (2). However, differences between the FNS 
and RNS solutions were observed when comparing the circumferential distortion as 
measured in terms of the AC«, engine face descriptor, Fig. (3). Both the FNS and RNS 
analysis indicated that the M2 129 inlet separated over the Mach number range from 0.1 
to 0.794. The location of separation, as determined from both the Full Navier-Stokes and 
Reduced Navier-Stokes solution techniques, is presented in Fig. (4), and suggests that 
some of the differences in the AC^ engine face distortion indicated in Fig. (3) can be at- 
tributed to differences in the prediction of flow separation location within the inlet duct. 
Figure (4) also indicates that the location of flow separation is effected by inlet throat 
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Mach number, and that no consistant statement can be made as to the prediction of 
separation using either FNS or RNS solution techniques. An examination of the pre- 
dicted separation location for both Full Navier-Stokes and Reduced Navier-Stokes ana- 
lyses relative to the experimentally measured separation point as well as a comparison of 
the experimentally measured inlet performance is presently underway. 

Presented in Figs. (5) through (8) are a comparison between the predicted engine 
face total pressure recovery map using both Reduced Navier-Stokes and Full Navier- 
Stokes solution techniques at the Test Case 3.1, 3.2, and 3.3 initial flow conditions. In 
general, there is strong similarity between the engine face flow field as computed by the 
two Navier-Stokes techniques, although the FNS solution at the Test Case 3.1 and 3.2 
flow conditions indicated a somewhat more developed recovery map. This was not true 
at the Test Case 3.3 initial flow conditions. The more developed engine face flow field 
computed by the FNS analysis is indicated in Fig. (4) as a higher DC^ engine face dis- 
tortion. At the Test Case 3.3 flow conditions, the underdevelopment of the engine face 
recovery map of the FNS solution, as compared to the RNS solution, is indicated as a 
lower DCfo distortion. 

The nature of the computed differences between the Full Navier-Stokes and Re- 
duced Navier-Stokes solutions technique is clearly indicated in Figs. (8) an (9), which 
presents the radial pressure ring distortion. Fig. (8), and the 60° sector circumferential 
ring distortion characteristics, Fig. (9). The radial pressure ring distortion as computed 
by both Navier-Stokes techniques is essential equivalent at the Test Case 3.2 initial con- 
ditions. This is also true at the Test Case 3.1 and 3.3 initial flow conditions. The primary 
differences in the prediction of flow separation as computed by Full Navier-Stokes and 
Reduce Navier-Stokes analyses shows as a difference in circumferential pressure dis- 
tortion as indicated in Fig. (9) and also Fig. (3). If agreement between solution tech- 
niques is measured in terms of standard engine face distortion descriptors, then the two 
Navier-Stokes analyses gave remarkably similar results for a inlet duct experiencing flow 
separation characterized as vortex lift-off. 

Presented in Figs. (10) through (13) are comparision between the static pressure 
and wall skin friction distribution along the 6 — 180 surface element for the Test Case 
3.1 initial conditions, Figs. (10) and (11), and at the Test Test Case 3.2 flow conditions, 
Figs. (12) and (13). The static pressures distributions as determined from the FNS and 
RNS solution technique are in surprising agreement. Fig. (10) and Fig. (12). The major 
difference between Full Navier-Stokes and Reduced Navier-Stokes solutions occure in the 
wall skin friction distribution upstream of flow separation. Figs. (11) and (13). i.e. the 
region of adverse pressure gradients in the first section of the S-duct. The skin friction in 
the separated region is characterized as very constant and very close to zero for both 
analysis techniques. Thus, the near wall velocity distributions within the region were 
vortex lift-off takes place are very small with very little variation relative to the velocities 
in the outer region of the flow. The differences in computed skin friction by the two 
Navier-Stokes solutions within the separated regions are very small. 

Figures (14) and (15) present the surface streamline pattern indicating the topology 
of vortex lift-off as given by the Reduced Navier-Stokes solution technique at the Test 
Case 3.3 initial condition. A very striking and significant feature captured by the analysis 
is the convergence of the limiting streamlines as an indication of three dimensional sepa- 
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ration taking place in this duct. Another important and striking feature is the symmetric 
pair of spiral nodes and pair of saddle points that were clearly captured by the 3D RNS 
analysis.. The topological patterns, as shown in the Figs. (14) and (15), also reveal the 
remarkable characteristic that the limiting streamlines forming the -spiral node enter only 
from downstream of the nodal point. The very familiar topological pattern is known to 
describe the important stage in the development of the pair of counter rotating vortices 
that form in the First section of turning resulting in vortex liftoff in the second section. 

CONCLUDING REMARKS 

A three-dimensional implicit Full Navier-Stokes (FNS) analysis and a 3D Reduced 
Navier Stokes (RNS) initial value space marching solution technique has been applied to 
a class of separated flow problems within a diffusing S-duct configuration characterized 
as vortex-liftoff. Both Full Navier-Stokes and Reduced Navier-Stokes solution tech- 
niques were able to capture to overall flow physics of vortex lift-off, however more con- 
sideration must be given to the development of turbulence models for the prediction of the 
locations of separation and reattachment. This accounts for some of the discrepancies in 
the prediction of the relevant inlet distortion descriptors, particularly circumferential dis- 
tortion. The 3D RNS solution technique using FLARE approximations adequately de- 
scribed the topological and topographical structure of flow separation associated with 
vortex liftoff. 
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Fig. (1) - Geometry definition for the M2129 intake duct. 


Mach 

Reynolds 

RNS Solution 

FNS Solution 

No. 

No. 

Xsep 

PR 

DC60 

Xsep 

PR 

ucou 

pm 

0.301xl0e6 

4.456 

WKM 


...... 

— 

— 


0.594x1 0e6 

4.897 



4.893 

0.996 

0.229 


0.873x1 0e6 

5.109 

WmMm 

0.254 

— — 

- " — -- 


1 0.412 

1.158xl0e6 

5.323 

0.990 

0.254 

4.456 

0.987 


wmm 

1.385xl0e6 

5.109 

0.985 

0.267 

— 



0.600 

1.577xl0e6 

4.675 

0.979 

0.297 

— 



0.700 

1.723xl0e6 

4.236 

0.971 

0.348 

— 

— - 


0.794 

1.848xl0e6 

3.794 

0.962 

0.416 

4.674 

0.958 

0.441 


Table I - Initial flow conditions and summary of inlet performance. 
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Total Pressure Recovery, Tt, t \PU 



Throat Mach Number, M-, 


Fig. (2) - Effect of inJet throat Mach number on engine face total 
pressure recovery for the M2 129 intake duct.. 



Throat Mach Number, Mi 


F jo (3) - Effect of inlet throat Mach number on DC W engine face 
° distortion for the M2129 intake duct. 
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Analysis 



Throat Mach Number, M-, 


Fig. (4) - Effect of inlet throat Mach number on the location of flow 
separation in the M2 129 intake duct. 



RNS ANALYSIS FNS ANALYSIS 

Fig. (5) - Engine face flow field for the M2129 intake duct. 

Test Case 3.1 initial conditions. 
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RNS ANALYSIS 


FNS ANALYSIS 


Fig. (6) - Engine face flow field for the M2 1 29 intake duct. 
Test Case 3.2 initial conditions. 
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Dimensionless Engine Face Radius, 6<f 


Fig. (8) - Radial pressure ring distortion characteristics. Test Case 3.2 
initial conditions. 



Dimensionless Engine Face Radius, 


Fig. (9) - 60° - sector circumferential pressure ring distortion characteristics. 
Test Case 3.2 initial conditions. 
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Wall Skin Friction Coefficient, Qx 10 2 « • Wall Static Pressure Ratio, pJPt ( 



Axial Centerline Distance, XjR; 


(10) - Static pressure distribution along the 0 = 180° surface element. 
Test Case 3.1 initial conditions. 



Axial Centerline Distance, X/R; 

Fig. (1 1) - Wall skin friction distribution along the 6 = 180° surface element, 
Test Case 3.1 initial conditions. 
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4 . 



Fig. (12) - Static pressure distribution along the 9 = 180° surface element, 
Test Case 3.2 initial conditions. 



Fig. (13) - Wall skin friction distribution along the 0 = 180° surface element, 
Test Case 3.2 initial conditions. 
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Fig. (14) - Surface streamlines pattern showing the topology of vortex liftoff. 
Reduced Navier-Stokes solution. Test Case 3.2 initial conditions. 



Fi«. (15) - Enlarged streamline pattern showing the topology of vortex liftoff, 
° Reduced Navier-Stokes solution. Test Case 3.2 initial conditions. 
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